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THE MODIFICATION AND APPLICATION OF KUO'S 
PARAMETERIZATION OF CUMULUS CONVECTION 
IN MIDDLE LATITUDES 

CHAPTER I 
INTRODUCTION 

Over many regions in the middle latitudes as well as in 
the tropics, cumulus convection plays an important role in 
vertical transport of heat, moisture and horizontal momentum, 
especially in the warm season (Palmen and Newton, 1969). In 
these regions , a large portion of the annual precipitation 
is derived from convective showers and thunderstorms; these 
storms also account for the greatest weather damage in many 
areas. Convective systems also are of particular interest 
because they manifest perhaps the most obvious interaction 
between disturbances of large and small (or medium) scales. 
However, because of insufficient resolution of the present 
observational network, this mechanism must be parameterized 
in terms of large-scale variables. In the past decade, many 
parameterizations of cumulus convection have been developed. 
Each scheme has served surprisingly well in the particular 
numerical model in which it has been employed. However, be- 
cause many investigators believe that cumulus convection is 
a dominant mechanism in the tropics (Riehl and Malkus 1958) , 
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most procedures have been developed for this region. We 
hardly can expect these schemes to fit our needs in middle 
latitudes without some modifications. 

Currently, two types of cumulus parameterizations are 
used. One, based on the hypothesis of convective adjustment 
(Manabe et , 1965) , is utilized in general circulation 
models and synoptic-scale prediction models. The idea of 
convective adjustment is based on the assumption that thermal 
convection develops when the lapse rate of temperature exceeds 
a certain neutral value. The other, based on the hypothesis 
of penetrative convection (Ooyama, 1964; Kuo, 1965; Arakawa, 
1969), is also utilized in general circulation models; but 
the scheme is primarily applied to the study of the develop- 
ment of tropical cyclones (Yamasaki, 1968; Ooyama, 1969). 
Penetrative convection is assumed to occur when cumulus clouds 
penetrate deeply into an unsaturated atmosphere in areas of 
low level mass convergence. Although both types of cumulus 
parameterization have practical significance, the latter 
appears to be ^more feasible for our purpose, because the 
hypothesis of penetrative convection is pertinent for middle- 
latitude deep convective systems and it can be used to simu- 
late cloud ensembles after some modifications are made. From 
the simulation of such an ensemble, one can obtain a great 
deal of information on the statistical properties of the 
cloud clusters and a deeper understanding of the dynamical 
mechanisms involved. 
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Based on the hypothesis of penetrative convection, the 
parameterization procedure consists of two significant mecha- 
nisms involved in the interaction between convection and the 
environment: 1) adiabatic warming due to the downward motion 

of the environmental air which compensates for the upward 
motion inside convective elements, 2) lateral mixing of cloud 
s\abstcince into the environment. Among such parameterization 
schemes, those developed by Kuo (1965, 1974) and Arakawa (1969, 
1974) are two of the most representative. Kuo's scheme is 
based on a non-steady deep cumulus model; the temperature 
environment, and the large-scale low-level convergence of 
moisture are the key indicators. Arakawa 's scheme is based 
on the concept of a balance of the vertical mass transport 
in the clouds and the environment. Although several authors, 
such as Ooyama (1971) and Fraedrich (1973, 1974), have de- 
veloped a theoretical framework to include these mechanisms 
and to discuss the relationship between these schemes, some 
confusion still exists (e.g., Ceselski, 1974). It has been 
stated that Kuo has neglected the heating by adiabatic com- 
pression of the slowly descending environment. Also, the 
reality of the horizontal mixing process involved in his 
scheme has been questioned. Recently, Kuo (1974) tried to 
clarify those misunderstandings; however, it is believed 
that if in Kuo's formulation the terms which account for 
departures from the area mean, as well as the precipitation 
terms , are examined carefully , a better explanation can be 
given. 
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Kuo's scheme is applicable in those situations where a 
deep conditionally unstable layer and large-scale low-level 
convergence are in existence. The former condition makes it 
possible for huge cumuli to penetrate into the upper tropo- 
sphere and the lower stratosphere, while the latter condition 
provides a lifting mechanism to trigger the convective insta- 
bility. Therefore, in some situations, although the computed 
fractional area of coverage may be near unity, only shallow 
or even no clouds may develop. Investigations of moisture 
convergence and its relationship to severe storm occurrence 
(e.g. , Sasaki 1973) indicate that unstable conditions, down- 
ward momentum transport, and other factors, which may be 
important to thunderstorm outbreaks, are implicitly shown 
in the moisture convergence patterns. Because we believe 
that large-scale moisture convergence in a conditionally 
unstable region is a key mechanism in the development of 
deep convection, we would like to examine the applicability 
of Kuo's scheme and to modify it so that the scheme becomes 
dynamically and practically plausible. 

In the following chapters, the AVE II data^ are utilized 
to perform the above purposes. In Chapter II, an objective 
analysis scheme is designed by a best combination of Inman's 
elliptical weighting function (personal communication) and 
Barnes' exponential weighting function (Barnes , 1973). In 

^Atmospheric Variability Experiments (AVE) conducted by 
NASA Marshall Space Flight Center, Huntsville, Alabama. 



Chapter III, a variational optimization scheme is developed 
to ensure accuracy of the divergence and vertical velocity 
fields. In Chapter IV, Kuo's parameterization scheme is 
examined and a modified parameterization scheme is proposed 
and investigated. It is believed that any plausible concepts, 
should be tested by a real situation. Therefore, in Chapter V, 
a case study is accomplished in order to check the variational 
optimization scheme, Kuo's parameterization scheme, and the 
modified parameterization of cumulus convection. 


CHAPTER II 


OBJECTIVE ANALYSIS 


In order to compute mass and moistxire convergence accu- 
rately, a variational optimization scheme has been designed 
to analyze the wind field. The scheme requires an objective 
analysis to interpolate values of wind to regularly spaced 
grid points from the irregularly spaced observation points. 
The weighting function utilized in the objective analysis 
should account for the observational resolution and real pat- 
terns of the data to avoid any unnecessary errors. Because 
the large-scale data field is anisotropic and nonhomogeneous , 
any isotropic or homogeneous weighting function will distort 
the data field and the unrealistic patterns may occur 
(McFarland, 1974). In this study, Inman's (1970) elliptic 
weighting function and Barnes' exponential weighting func- 
tion (Barnes, 1973) are combined to produce 


W(r,k) = exp {- — - Ijcos^e 




r 




( 1 ) 


where 


k2 = - 


In W, ' 

D 


r is the distance between the grid point and the observation 
point, b is the length of the minor axis of an ellipse of 
constant weight, 0 is the angle between the wind direction 
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and the direction from the grid point to the observation point, 
«b is the weight when r is equal to b and 0 is 90 deg, t is 
the time difference from a reference time, v is a constant and 
a is the ratio of the length of major axis to that of the minor 
axis b of an ellipse of constant weight* In this study, a is 
chosen as 

ot = a ^ , (2) 

'^max 

where V is the wind speed at the observation point, V is 
the maximum value of the observed wind, and a is a specified 
constant. 

After the computation of the weighting function from (1) 
at each grid point, the value of a quantity x^j the grid 
point (i, j) is assigned as 






(3) 


where Xjjj is the value of the quantity x the m-th observa- 
tion point. 

In order to better describe the meteorological fields 
near a frontal zone, the interpolation scheme is utilized 
twice. In the vicinity of a front, the stations whose winds 
make a large angle with the wind at a grid point may be separated 
from the grid point by the front. After the first conventional 
interpolation, the wind is assigned at each grid point and the 
angle between the wind at a station and the wind at a grid point 
can be determined. During the second interpolation, any 
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stations whose wind make an angle with, the wind at a grid 
point greater than 90 degrees are not utilized in determining 
the analysis at the grid point. This procedure tends to 
conserve discontinuities, and is valuable especially near 
frontal zones. 
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CHAPTER III 


u 


VARIATIONAL OPTIMIZATION OF WIND FIELD 

In a region of severe convection the vertical velocity 
is not necessarily zero at an arbitrary level near the cloud 
tops. However, over a very large region the average vertical 

3^s 

velocity should be zero, if is negligible, to insure con- 

servation of mass. The integral constraint which requires the 
mean vertical velocity at the top over a large area to be a 
predetermined value (usually zero) is called the global boun- 
dary condition. In this chapter, Sasaki's variational opti- 
mization approach is utilized to optimize the entire wind 
field so that the integrated continuity equation and global 
boundary condition are satisfied. 

In the variational formulation, the upper boundary con- 
dition on the vertical velocity at the top of the convective 
layer should be specified in order to solve the associated 
Euler-Langrange equations. However, this condition is unknown. 
McGinley (1973) assumed that it was zero everywhere. O'Brien 
(1970) corrected the vertical velocity at the top, which was 
obtained from the continuity equation and the observed wind 
field, by the use of the global boundary condition as a strong 
constraint. However, because the areal mean of the vertical 
velocity at this level, is. usually small, the correction is 
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quite small at most grid points. Fxarthermore , since the cor- 
rection is a constant over a horizontal plane (or constant 
pressure surface) , it does not change the vertical velocity 
patterns. This is desirable only when the vertical velocity 
patterns obtained from the continuity equation and the observed 
winds are highly reliable. Indeed, McGinley's top boundary 
condition is an extreme case of O'Brien's, and is too arbitrary, 
especially in the vicinity of severe convective systems. 

The errors which appear in the vertical velocity field 
obtained from the continuity equation and the observed wind 
can be classified as "systematic errors" and "random errors". 

The systematic errors occur primarily because of incon- 
sistencies between the observed fields and the dynamical 
model considered. The random errors may be introduced because 
of inaccuracies of measurements , spatial irregularity of 
observation points , and by interpolation of values from sta- 
tions to grid points. Sasaki's variational optimization 
scheme provides an important advantage by incorporating 
dynamic, kinematic, statistical and other conditions in data 
management. Although the random errors can be suppressed by 
the use of varied filters, the filtering should be incorporated 
simultaneously in the variational formulation to ensure that 
the filtered fields satisfy all constraints. 

In the following paragraphs, a variational optimization 
scheme is developed such that systematic errors are suppressed 
by the use of two strong constraints, i.e., the integrated 
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continuity equation and the global boundary condition; the 
random errors are filtered by the use of a low-pass filter. 
The functional is defined as 


I = ///{(u-u)^+ (v-v)^ + + 3^ (Vo)^* Voi^) ^}dp dy dx. 

( 4 ) 


X y p 


Here, two strong constraints are utilized: 

,P4 


«t = ^s - / dp , 


( 5 ) 


// 


co^ dx dy = 0 , 


( 6 ) 


where u, u, v, and v are observed and optimized wind compon- 
ents along the x and y directions, respectively, co^ is the 
optimized p velocity at the top of the domain, a and 3 are 
specified parameters, p and p. are pressures at the surface 
and at the top of atmosphere, respectively, and 


/Pt . 

u. = m - / V*V dp , 


oj = p 

s s ^ 


s ' 


( 7 ) 


( 8 ) 


where ^ is the observed surface wind and V is a two dimen- 
s 


sional Laplacian differential operator, 

y = ^ ) 

' 3y^ * 


( 9 ) 


After taking the variation of the functional I and 
setting it to zero , the associated Euler- Lagrange equations 
are obtained. These may be written as 
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3A 


O ./V 

u = u + 3 ^ , 

(10) 

3A, 

V = V + . 

(11) 

<jy 

^1 -2- 
“t p -p^ ^2 - a , 

(12) 


where the associated boundary conditions are 

> 1=0 , 


(13) 


and 


•n»Vw^ = 0 / (14) 

where n is an unit normal vector along the boundary. The 
derivation of the above equations and the associated boundary 
conditions are discussed fully in Appendix C. 

In order to better understand the physical meaning of 
this formulation and the associated boundary conditions, let 
us define the correction velocity, ^* , as 

^ ^ , (15) 

where ^ = (u,v) , ^ = (u,v), and = (u*,v*). Substitution 
of (15) into (4) , (10) , and (11) leads to 

X =///<. *‘V* + + 3^(Vw^*Vw^)} dp dy dx / (16) 

X y p 

and 

. (17) 


Eq. (17) shows that the Lagrange multiplier is the correc- 
tion velocity potential; i.e., the correction velocity field 


is irrotational. The first two terms in (4) are the weighted 
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correction kinetic energy of the region considered, and a is 
the weighting factor. The third term in (4) is a weak con- 
straint, i.e. , 

- 0 ^ (18) 


Indeed, |Vw^| is a part of the horizontal vorticity associated 
with 0 )^. It reflects the spatial distribution of areas of 
upward and downward motion. In terms of the vertical correc- 
tion velocity lo*, (18) can be written as 

|Vw* + Vai^l = 0 , (19) 


or 


3w* 

<55T * 


3w 


3x 


"^)2 + 


3w. ^ 

( — £ + — ^ 
^3y 3y ^ 


0 . 


( 20 ) 


Since the horizontal gradient of the vertical velocity is 

zero only when 


3w* 

3x ~ 3x ' 


and 


3ujJ 3u, 

^ " 3y“ ' 


( 21 ) 


( 22 ) 


(19) indicates that (18) is satisfied by introducing a hori- 
zontal vorticity |Vw*| opposite to the original horizontal 
vorticity associated with Vw^. Since the horizontal vorticity 
indicates the spatial distribution of upward and downward 
motion, any reduction of the horizontal vorticity suppresses 
the vertical motion. In other words, this procedure smooths 
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the vertical motion field. The parameter 3 is selected to 
accomplish this purpose. In this study, the value of 3 is 
assigned to suppress the amplitude of waves whose wave length 
is less than the distance between adjacent maxima in the 
large-scale to field. 

The boundary condition (13) implies 

^ = 0 ^ (23) 

along, the boundary. Indeed, it is the no-slip boundary con- 
dition for the correction velocity field. This condition 
ensures that the optimization does not change the circulation 
of the region considered and is consistent with an irrotational 
correction velocity field. Furthermore, since co,^. is unknown, 
the boundary condition (14) permits to be free on the 
boundary to allow a wide range of possible patterns. 

Eqs. (5), (6), (10), (11), and (12) are five equations 

for the five unknown, u, v, and Their solutions, 

derivations and associated response functions are discussed 
fully in Appendix C. Also, in Appendix D, it is shown that 
O'Brien's and McGinley's formulations are two special cases 
of this formulation. 

A case study, utilizing NASA AVE II data, is discussed 
in Chapter V. The results are very encouraging. 



CHAPTER IV 

MODIFICATIONS OF KUO'S PARAMETERIZATION 
OF CUMULUS CONVECTION 

A Procedure for Evaluating Cumulus 
Parameterization Schemes 

Described in the following paragraphs is a procedure 
for examining and evaluating the mechanisms involved in Kuo's 
(1965) parameterization scheme and in the modified scheme to 
be specified later. The method is based on the concept that 
utilization of an accurate expression for the latent heat 
released by cumulus should lead to good evaluations of other 
characteristics of the cloud ensemble. It will be shown below 
that when Q^, the latent heat released by cumulus, is expressed 
in terms of certain cloud ensemble properties, the cloud en- 
semble properties may be determined from the large-scale 
budget equations if the time derivatives of dry static energy 
and mixing ratio in the environment is specified. In other 
words the cloud ensemble properties computed from this proce- 
dure should directly reflect the quality of the expression, 
for the latent heat released, which is utilized in the para- 
meterization scheme. Furthermore, if the large-scale budgets 
and the parameterization scheme are both, perfect, the latent 
heat released, as computed from the parameterization scheme 
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and as determined using the evaluation procedure, should be 
exactly the same. In the following discussion, any quantity, 
except the vertical velocity, in the environment is approxi- 
mated by the area-mean value of that quantity (for more 
detail, see Appendix A). 

First, let us define the dry static energy s and the 
m.oist static energy h as 


s = c T + gz , and 
P 

h = CpT + gz + Lq , 

where T is temperature, z is height, q is mixing ratio, Cp 
is the specific heat of dry air, and g is the acceleration 
of gravity. From (A. 2) and (A. 3), the heat and moisture 
equations for the large-scale motion can be written as 


0i - Qr ^ al - °R = «o 


- Qj = L g = + L - 5)1 , 


( 25 ) 


where s , q , s and q are dry static energies and mixing 
o o 

ratios in the cloud and in the environment, respectively, 

is the heating rate by radiation, L is the latent heat of 

condensation, p is pressure, Q is the latent heat released 

c 

by subgrid scale convection and M is the cloud mass flux de- 

o 

defined as = -am^, where a is the fractional area covered 
by the clouds, and m is the pressure velocity inside the clouds. 

w 

Also, the saturation moist static energy of the environ- 
ment f h*, is _ _ _ _ 

h* i CpT + gz + Lq* = h* , 
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where q* is the saturation mixing ratio at temperature T. 
Then, after Arakawa (1969) , we have 

® - h*' ' 

(h^ - h*) , (27) 

where y = ^(^^)_ • 
p 3T P 

If is expressed in terms of the other variables, 

combination of (24) - (27) gives four equations for the four 

unknowns, s , q , M , and h . They can be solved provided 
o c c o 

that the observed large-scale heat and moisture budgets 
over the area considered are known. The associated boundary 
conditions are 



- °R> = - -r'^cb - =b> ' 


- Qo = LE ^ L(q . - q, ) , 

g 2 s g ^cb ^b ' 


V^S = - i^)] , 

where the subscript b denotes values at cloud base and T , 

s 

To, q^, and q^ are the temperatures and mixing ratios at the 
surface and at a level above which is representative of the sur- 
face boundary layer. Let us define Y as - w'h' . Then, adding 
(24) and (25) and integrating the resultant equation, we get 

Y(P) = - Qj,) dp , 

= Me (he - H) , if p < pj^ , (28) 


I 
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where is the pressure at the top of atmosphere (100 mb in 
this study); Y(p^) is assumed to be zero. Utilizing (28), 
the boundary conditions are explicitly expressed as 


where 


and 


M V. = 
cb 


_ ®2 


h* - h. 
b b 


cb 


^b M 


cb 


= g(s^ + LE^) - Ap(Q^ - Qjj - Qj) , 


where p 
^s 


S + LE = Y(p ) , 

s s ^^s ' 

is surface pressure. 


(29) 


(30) 


An Examination of Kuo's Expression 
for the Release of Latent Heat by Ciomulus 


In this section, following Kuo (1965) , Q is specified 

c ^ 


as 

% = B(s^ - i) , 

which is (A. 13) expressed in terms of variables defined in 
the preceding section. The cloud base is assumed to be at 
the lifted condensation level representative of the surface 
layer, and the top of cloud is assumed to be at a level where 
the temperature of the clouds is equal to that of the envir- 


onment . 
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B is expressed according to (32), where b, and 


are defined by (A. 12), (A. 9), and (A. 10), respectively. Sub- 

stitution of (A. 12) and (A. 10) into the expression for B 


leads to 


B 


“i /( 


Q, ^ 
6 2 g 


Because b is computed from the local time change of q, from 
the parameterization point of view b is unknown and must be 
specified experimentally. In order to avoid additional error 
introduced from the specification of values of b , in this 
study b is computed from 9q/3t according to (A. 12). The 
same values of b are used in the computation of the release 
of latent heat in the parameterization schemes. 

After elimination of s , q and h from (24) - (27) , 

c c c 

an ordinary nonlinear differential equation is obtained; 


dp 


F - G/U , 


(31) 


where 


U = M (fc^) 
c'l+Y ' 


t 


F 



B(h-h*) 3 . Y , 

1+Y 3p^l+Y^ ' 


G = §I_(h-h*) , and 

(1+Y)^ 

B = j^(l - b) , 

“l 


(32) 


where Y is determined from (28) as a function of Q^, Q^, and 
Qj^, and b is defined by (A. 12). 

Eq. (31) can be reduced to Abel's equation of the first 
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kind and can be solved analytically if the ratio of the two 

coefficients, G and F, is constant, or it can be solved 

numerically even if the coefficients are variable. is 

determined from solutions of (31). Then, according to (28) 

h can be determined by 
o 


^ + h . 

C M 

c 

Next s^ and are computed utilizing (26) and (27) , respec- 
tively, and is determined by 

®c " ® - 

T = + T 

C C_ ^ ^ • 


Finally, may be calculated from (A. 13). 

Employing the procedures described in the preceding 

paragraphs, calculations of M , h , T , q , Q , and heights 

C 0* o 

of cloud tops have been made utilizing NASA AVE II data. 
Results of these calculations are shown in Figs. 33-36. Al- 
though radar charts prepared by the National Weather Service 
(2035 GMT 11 May 1974) show that the cloud tops were at about 
275 mb in the vicinity of the test area, the computed cloud 
tops are at 556 mb. M computed from (31) is unrealistically 
large. For example, the values are about an order of magni- 
tude larger than those calculated by Lewis (1975) for a 
thunderstorm system in central Oklahoma. Also, the latent 
heat released is about an order of magnitude smaller than 
required by the large-scale budget equations. All results 
show that the properties of the cloud ensemble, computed on 
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the basis of Kuo's expression for the latent heat released, 
are questionable. 

In his tropical prediction model, Krishnamurti (1969) 
treated the terms which account for departures from the area 
mean (i.e. vertical diffusion terms) and the precipitation 
terms, separately. The precipitation terms were computed 
utilizing an extension of Kuo's scheme, and the vertical 
diffusion terms were computed utilizing eddy exchange co- 
efficients; the results seem satisfactory. However, other 
problems still exist. Reed and Recker (1971) found evidence, 
in a composite Pacific tropical wave, of a heating maximum 
centered near 400 mb. As mentioned by Ceselski (1974), Kuo's 
scheme will yield an upper tropospheric heating maximum only 
if unrealistically deep and hot clouds are assumed. Clouds 
that include entrainment and extend only to 200 mb or below 
will generally have maximum (T - T) in the lower troposphere 
(Ceselski, 1973). This is particularly true if the distur- 
bance in question is relatively warm in the upper troposphere, 
as is often the case. 

Recently, Edmon and Vincent (1975) have employed Kuo's 
scheme, as well as the modified scheme by Krishnamurti, to 
calculate the latent heat released in a case involving the 
interaction of tropical storm Candy with an intensifying 
extratropical baroclinic frontal system (0000 - 1200 GMT 
25 June 1968). Comparisons between the convective latent 
heat released and actual precipitation rates show, in general. 
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that Kuo's scheme and Krishnamurti ' s scheme fail. They state 
that there was very little convective latent heat released 
in the computational region, even though there were widespread 
thunderstorms during the analysis periods. Therefore, it is 
necessary to modify Kuo's formulation so that the modified 
scheme may be applied in middle latitudes. 

Modified Parameterization of Cumulus convection 
Before discussing any modifications, a few statements 
about the cloud ensemble model proposed in this study are appro- 
priate. Here we assume that the hydro- thermodynamic fluid 
properties in area a are the same as those of the environment 
before the clouds form. Later, clouds form and the whole area 
a is covered by clouds, and the hydro- thermodynamic fluid 
properties in a are changed from values typical of the envir- 
onment. During the formation period, entrainment and detrain- 
ment, condensation, evaporation, and induced subsidence in 
the environment may occur. The formation period can be 
divided into two stages; in the first stage moisture is 
supplied to the cloud and in the second stage condensation 
occurs. The conservation law for water vapor implies that 
there are two sources of water vapor available to produce 
clouds in a layer. One is the large-scale moisture supply 
(i.e., the apparent moisture sink) through the lateral 
boundaries, and the other is vertical moisture convergence 
inside the cloud across the upper and lower bo\indaries of 
the layer. In other words, during the first stage, for the 



modified scheme, the rate of increase of mixing ratio in a 
layer in the cloudy region is formulated as 


B = 

q 






(33) 


The second term on the right-hand-side of (33) allows the low 
level moisture convergence to produce a high level release 
of latent heat. If (33) is integrated from the top of atmos- 
phere down to the surface, we have 


I f \ 

g^Pt q 


= --f 

g 


Pc; 

3t 


dp - 




qV dp + E 


= (1 - b)Mj^. 

The above expression is the same as that used by Kuo to deter- 
mine the rate at which moisture is supplied to form clouds. 

During the condensation stage, a part of the moisture 

supplied to the cloudy region is utilized to increase the 

mixing ratio from q to q^, and the other part is condensed 

and latent heat is released. Since the moisture required to 

produce a unit mass of cloud air (i.e., to increase the 

mixing ratio and temperature from q ,T to q ,T ) in a layer of 

o o 

thickness, (5p is ^(h^ - ii) , for the modified scheme, the pro- 
duction rate of cloudy air in a layer is expressed as 


B - 


B. 


(h^-h)/L 


if B >0 

q 


0 


if Bg J< 0 . 
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On the other hand, the conservation law of energy implies 
that the latent heat released in a layer can be divided into 
two parts. One part is utilized to heat the air inside the 
cloud in that layer from T to and the other part is trans- 
ported to adjacent layers by the vertical motion in the 
cloudy region. Therefore, for the modified scheme, the 
latent heat released is formulated as 

= B(s^ - i) - , (34) 

where - 'the vertical divergence of dry static 

energy in the cloudy region. The latent heat released is 

zero if Q is negative; the level where is zero is assumed 
c c 

to be the top of clouds. In other words, the latent heat 

released is assiomed to be zero at the cloud top. If M in 

c 

(34) is zero,. and B is replaced by (1 - b)Mj^/M^, (34) is 

reduced to Kuo's expression for the latent heat released. As 
shown by Fraedrich (1973) , Kuo did not account for the verti- 
cal mass flux inside the cloud. Comparison of (34) with (A. 13) 
indicates that a significant modification is that the effect 
of cloud mass flux on the latent heat released is included 
in the modified formulation. 

A comparison of the modified scheme with the cloud 
ensemble model proposed by Yanai et al. (1973) , may clarify 
the physical processes allowed in the modified scheme. Accord- 
ing to Yanai (1973) , if the cloud ensemble is assxxmed to be 
in quasi-equilibrium with the large-scale forcing mechanism 



25 


(Arakawa and Schubert, 1974) , the budget equation for heat 
inside the clouds may be written as 


fp<Mc=c> IH^ - f i=oi ^ = 


(35) 


where e . , D . and s . are entrainment , detrainment and dry 

l' X Cl 

static energy of the i-th cloud, respectively, and C is the 
condensation inside the clouds. Since 


= L(C-e) , (36) 

where e is the rate of evaporation of liquid water detrained 
from the clouds, combination of (34), (35) and (36) implies 


B(s -s) = - Ze.s + ZD.s . - Le . (37) 

c ^ X 1 Cl 

Eq. (37) shows that entrainment, detrainment, and evaporation 
are parameterized in terms of the temperature difference be- 
tween the cloud and the environment and the production rate 
of cloud air. Furthermore, substitution of (34) into (24) 
leads to 


_ 3M 

Ql - Qr = - “c If - ^ 8?^ ^ ®(^c - • 


The above equation shows that the apparent heat source for 

the large-scale motion field, in addition to Qj^, consists of 

cui adiabatic warming due to downward motion outside the 

clouds which compensates upward motion (positive M ) inside 

c 

the clouds , a cooling due to the vertical mass divergence 
outside the clouds which compensates the vertical mass con- 
vergence inside the clouds, and a warming directly due to 
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the increasing dry static energy inside the clouds. 

For the modified parameterization scheme, in order to 
determine the latent heat released utilizing (34) , the cloud 
mass flux has to be determined first. A simple two- layer 
model, which consists of a primary and a detrainment layer, 
is suggested to accomplish this purpose. Except near the 
cloud top where a strong detrainment layer exists due to the 
strong negative buoyancy above the cloud top, the vertical 
mass convergence inside the cloud in the primary layer is 
assumed to be proportional to the mass convergence of the 
large-scale motion and to the fractional increment of mois- 
ture inside the clouds; i.e.. 


^ M = XM^ , where 
3p c c 

q -q 

X = - + K2(-|— ) . 


(38) 


(39) 


Indeed, (38) is similar to the conventional entrainment model. 

In the detrainment layer near the top of the cloud, the 
cloud ensemble model is described by the following three 
equations ; 

0 = - D + , (40) 


0 


Ds + ^(s M ) + LC 
c 3p c c 


(41) 


0 = - Dq + (q M ) - C . 
^c 3p ^c c 

The latent heat released is specified as 


®c = - = 


(42) 


LC . 


(43) 
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Eqs. (40) , (41) , and (42) are similar to those of the cloud 
ensemble model proposed by Yanai et (1973) except that 

here entrainment is assumed to be negligible compared to 
detrainment. Eq. (43) shows that production of cloud air is 
neglected. Evaporation also is neglected because the liquid 
water content of the clouds is very small at high levels. 

The solutions of (40), (41), and (42) are 


3M 
c 

3p 

LC 


3p c 


= D , 


■ ^c 3p ®c ' 


= 0 . 


(44) 

(45) 

(46) 


Eqs. (44), (45), and (46) are three equations for four 

unknowns, T , M , C, and D. In this section, one more equa- 
o o 

tion is introduced by assuming that the detrainment, D, is 
constant in the layer such that the cloud mass flux is 
zero at the top of the clouds; i.e.. 



where Ap is the depth of the detrainment layer, and M* is 

o 

the cloud mass flux at the bottom of the detrainment layer. 

In order to make sure that the modified parameterization 
scheme is reasonable, the procedures described in the first 
section of this chapter will be utilized; results of this 
investigation will be discussed in Chapter V. Also, the 
modified scheme will be compared to Kuo's (1965) scheme. 



CEAPTER V 


TESTS 

Optimization of the Wind Field 
All data are from the second NASA Atmospheric Variabi- 
lity Experiment (AVE II) except for the radiation data. There 
were fifty-four rawinsonde stations participating in the AVE II 
Pilot Experiment as shown in Appendix B. Soundings were made 
at three-hour intervals at each station beginning at 1200 GMT 
11 May 1974, and ending at 1200 GMT 12 May 1974. The data 
were obtained during a period when convective activity was 
present, large horizontal temperature gradients were evident 
and rapid changes in weather patterns were occurring. The 
data area is over the eastern United States east of approxi- 
mately 105° W longitude. Radar data were obtained from ele- 
ven stations located near the center of the observational 
area, and as much data as possible were collected from the 
NIMBUS 5, NOAA-3, ATS-3, and DMSP (DAPP) satellites. 

The radiation climatological data are from Rodger's 
results (1967) , which are available in the region 0-70 N, 
1000-10 mb, and for the months of January, April, July and 
October. Because the AVE II pilot experiment was held in 
May, an average of the data for April and July will be 
utilized in this research. 


28 




The synoptic situation of 2100 GMT 11 May 1974 is shown 
in Figs. 2-8. We choose this particular time because a deep 
conditionally unstable layer exists over a wide area of the 
southern United States. There is a cold front across the 
central part of the coxintry and a warm front extends through 
the northern states. Three lows and two precipitation areas 
appear on the map at the analysis time. Also, there is a 
deep trough in the upper levels extending from north to south 
across the central United States. 

The grid system is shown in Fig. 1; the polar stereo- 
graphic projection is utilized. The standard latitude and 
the standard longitude are 60“N and 100®W, respectively. The 
11 by 13 computational grid is oriented so that the y-axis is 
perpendicular to the standard longitude. The map scale is 
1:15,000,000, the upper- left-hand corner grid point is located 
at X = 12.6 in and y = 0.97 in (the origin is at the North 
pole), and the grid interval is 190.5 km on the image plane. 

The interpolation scheme is described in Chapter II. 

The parameters a, k and v in (2.2) and (2.1) are / 8, 0.9772, 
and 2, respectively. During the second interpolation, the 
weighting function is assigned as W 2 = cos 4», where tp 
is the angle between the wind at observation stations and 
the wind at grid points, and is the weighting function of 
the first interpolation determined from (2.1). After in- 
terpolation of each variable (i.e. , temperature, geopotential 
height, wind velocity, and mixing ratio) to grid points, mass 
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divergence and moisture divergence are computed utilizing a 
fourth-order finite-difference approximation. The vertical 
p velocity f w, is determined at each grid point, by integration 
of the continuity equation from the surface; its surface 
value is determined from (8). Results are shown in Figs. 9-17. 
The values of tu at the surface are very small (10 ^ mb sec ; 
the patterns reflect the synoptic situation and orographic 
lifting. Because of the accumulative errors in the diver- 
gence fields, the o) fields at higher levels are highly un*' 
realistic, both in pattern and in magnitude. For example, the 
01 patterns don ' t clearly reveal the areas of upward and down- 
ward motion associated with the two troughs in the upper levels, 
and the magnitudes of the oi values are incredibly large, 
especially at the 100-mb level. 

The optimization of the wind field is performed utilizing 
0.1 and 0.3, respectively, for the values of the parameters a 
and 3 in (4). Results are shown in Figs. 18 - 25. From 
600 mb up to the 100-mb level, the optimized oi fields clearly 
reveal a major trough across the central United States and a 
minor trough across the southwestern portion of the country. 

To the east of the major trough, strong upward motion exists 
and to the west of the major trough strong downward motion is 
evident in the optimized oj fields. Also, there is a small 
area of upward motion associated with the minor trough in the 
southwestern United States. All these features were ambigu- 
ous before the optimization was accomplished. Also, the 
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values of lu are reduced to reasonable magnitudes by the opti- 
mization. Compared to the National Weather Service radar 
chart for 2035 GMT 11 May 1974, the area of upward motion 
coincides with the region of deep convection. 

Tests of the P a raraeterization of Cumulus Convection 
In the following paragraphs, Kuo's scheme and the modi- 
fied scheme are examined by utilizing the procedures discussed 
in Chapter IV. Eqs. (24) - (27) along with the associated 
boundary conditions are used in this examination. Since 

these four equations contain five unknowns, s , q , h , M 

c c o c 

and Q , an additional equation must be obtained from the for- 
o 

mulation of latent heat released. Different formulations of 
provide different solutions. When Kuo's scheme is examin- 
ed, is formulated according to (A. 13), and when the modi- 
fied scheme is investigated, is formulated according to 
(34). The solutions include the latent heat released, the 
height of the cloud tops , and four bulk properties of the 

cloud ensemble, s , q , h and M . These bulk properties, 

o o o o 

the latent heat released, and the height of the cloud tops 
reveal physical insights of the mechanisms involved in each 
parameterization scheme. 

After the evaluation procedure is completed, the latent 
heat released by cumuli will be determined directly, utiliz- 
ing both parameterization schemes. The computed latent heat 
released will be compared to the observations, i.e.. 


to the 
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precipitation which is required by the large-scale budgets, 
and to the solutions obtained by utilizing the evaluation 
procedures specified in Chapter IV. 

The tests are performed at grid point (8,9) where a 
deep unstable layer exists from the lifted condensation level 
of the surface layer up to the 100-mb level. The horizontal 
diffusion terms V*V's' and V»^'q' in (4.1) and (4.2) are com- 
puted by the use of nonlinear horizontal diffusion coeffici- 
ents similar to those used by Smagorinsky e_t a)^. (1965) and 
Manabe et al. (1965); i.e.. 


VV's' = V* (K 7s) , 
s 


V»^'q' = V* (K- Vq) , 

"q = "s = I ' 


where 


9u _ ^ 
3x 3y ' 


9v . 9u 
8x 3y 


is the Karman constant (0.4 in this study), and A is the 
grid length. The results show that values of these horizon- 
tal diffusion terms are one or two orders of magnitude smaller 
than those of the apparent heat source, Q^, and the apparent 
moisture sink, Q^. The local time changes of the dry static 
energy and the mixing ratio are determined from observations 
at 1800 GMT and at 2400 GMT 11 May 1974. Results show that 
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these two terms are not negligible. The divergence field 
cind the u field are optimized according to the procedures 
described in Chapter III and in an earlier section of this 
chapter, and are shown in Figs. 26 and 27. The lifted con- 
densation level is computed according to Inman's approximate 
formula (1969). The temperature T, the mixing ratio q, the 
dry static energy s, the moist static energy h, the satura- 
tion moist static energy h*, the apparent heat source 
the apparent moisture sink Q^f the radiation and the 

derived vertical eddy heat flux Y are shown in Figs. 28-32. 
Because of the large horizontal gradients of moisture and 
temperature in the middle latitudes , the apparent heat source , 
the apparent moisture sink, and the derived vertical eddy 
heat flux are about three to five times larger than those 
of the tropics as determined by Yanai et aJ^. (1973). 

In order to obtain a better comparison between the ob- 
servations and the results computed from the parameterization 
schemes, the required precipitation is determined from the 
large-scale heat and moisture budgets as follows. Let us 
assvune that the vertical eddy heat flux and the vertical 
eddy moisture flux is zero at the cloud tops, i.e.. 


M (s - s) = 0 , 

P=Pt 


Mc(q^ - q) 


P=P, 


0 
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Then, integration of (24) leads to 
1 

I =s ' 

where is the total precipitation inside the clouds , and 
Sg is the rate of transport of sensible heat from the surface. 
Since the observed tops of clouds are near 275 mb, computa- 
tions based on the large-scale heat and moisture budgets 
indicate that the required is 0.08 in hr This amount 

of precipation is reasonable compared with observations from 
the National Weather Service radar charts (2035 GMT 11 May 
1974) and Service A teletype reports. 

In the evaluation of Kuo's parameterization scheme 
is determined from (31) and B is specified by (32) . Because 
the production rate of cloud air, B, depends on M^, the 
equations are nonlinear integral-differential equations. A 
Runge-Kutta method and an iterative technique are utilized to 
solve the set of equations. The solutions are shown in 
Figs. 33-36. From these results and (47) the precipitation 
rate is 0.00966 in hr and the cloud top is at 556 mb. 
Compared to observations , the computed precipitation is one 
order of magnitude smaller than that observed, and the top 
of the cloud is about 300 mb lower. Furthermore, the cloud 

mass flux M and the residual mass flvix in the environment, 
o 

M, are unrealistically large. These unsatisfactory results 
may be understood by the following discussion. 



Let us rewrite (24) and (25) as follows: 

°o = °1 - °R - 

Q^A = Qj/L + - 5)1 . (49) 

The first two terms on the right-hand-side of (48) indicate 
that a part of the latent heat released in a layer is utilized 
to heat the atmosphere in that layer, and the third term de- 
notes that the other part of the latent heat released is 
utilized to heat the air in the adjacent layers by means of 
the vertical motion inside the clouds. The first term on the 
right-hand-side of (49) indicates that a part of the precipi- 
tation in a layer comes from condensation of water vapor in 
that layer, and the second term shows that the other part of 
the precipitation comes from the moisture of adjacent layers 
by means of the vertical motion in the clouds. In other 
words, latent heat is released as condensation occurs and the 
dry static energy in a layer is increased. At the same time 
the vertical eddy flux terms account for the redistribution 
of moisture and dry static energy. However, because of large 
horizontal gradients of temperature and moisture in the middle 
latitudes, the apparent heat source, Q^, and the apparent 
moisture sink, Q„, are very large and (s - s) and (q - q) 

A w C 

are very small. The vertical eddy flux terms may transport 
the required energy and moisture only when the vertical mass 
divergence (or convergence) inside the cloud x— ^ 


/ 


is extremely 
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strong. In other words , the cloud mass flux must increase 

rapidly with height. According to (26) and (28), the cloud 

temperature can be expressed in terms of M , Y, and large- 

c 

scale variables as 

’’c = ^ + Cp(U) 'i; - ^<5* - 5)] • 

From the above equation, it can be seen that as M increases 
the cloud temperature, decreases. When is sufficiently 
large that 

“c ^ ' 

the cloud temperature will be less than that of the environ- 
ment. In other words, the rapid increase of M with height 
decreases the cloud temperature and, therefore, suppresses 
the vertical development of the cloud. This phenomena is 
verified by the low tops of clouds and the large cloud mass 
flux in the solutions of the examination of Kuo's scheme. 

The discussions indicate that a suitable modified scheme 
should consist of a better mechanism to transport energy 
and moisture vertically. 

The equations and associated boundary conditions util- 
ized in testing the modified parameterization scheme are the 
same as those used in testing Kuo's scheme except that the 
latent heat released is formulated according to (34). The 
equations are nonlinear ordinary differential equations, and 
a Runge-Kutta method is utilized to solve them. The solutions 
are shown in Figs. 37-39. From these calculations and (47) , 
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the precipitation rate is 0.075 in hr the tops of clouds 
are at 260 mb, and the maximum latent heat released is at 
450 .mb. These results agree very well with the observed 
values. 

In order to better explain the improvement of the modi- 
fied scheme, let us decompose the total latent heat released, 
as given by (34) , into three components as shown below; 


Qc = °2 




(h^-h) h^-h 


L(s -s) j, - 

^ ^ - («c‘3cJ - Ip <^c^c> • 




Each term on the right-hand- side of the above equation 
is evaluated and the results are plotted in Fig. 41. The 
first component is due to the large-scale moisture conver- 
gence, the second is due to the vertical transport of mois- 
ture, inside the cloud, and the third represents the vertical 
divergence of the dry static energy. These three components 
are nearly equal in magnitude. Therefore, the latent heat 
released in a layer may not be confined to that layer; the 
vertical transport of released latent heat is a significant 
mechanism involved in cimulus convection. Because of the 
third term, the maximum of the latent heat released is at 
the 450-mb level. In order to gain insight of the physical 
meaning of this component, let us consider a special case. 

If the apparent heat source, is negative such that the 

rate of increase of moisture in a layer inside the cloud, B , 
is zero, i.e.. 


9 


(50) 
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and (34) becomes 


°o = - lp<Vo' ' 


( 51 ) 


In this situation, although the production rate of cloud air 
is zero, latent heat is still being released because of the 
nonzero vertical mass flux inside the clouds. The cloud acts 
as a machine which pumps the moisture and the dry static 
energy from the lower levels up to higher levels or vice 
versa. Also, because the freezing process may occur at the 
upper levels, the maximum latent heat released may be at the 
higher levels. 

So far, examination of the modified scheme by the use 

of the large-scale budget equations shows that the modified 

formulation of latent heat released can better simulate the 

physical mechanisms involved in the release of latent heat. 

However, from the parameterization point of view, the cloud 

mass flux, M , and the cloud temperature have to be deter- 
o 

mined first in order to utilize the modified scheme to 

compute the release of latent heat. In this section, is 

computed according to (38) and (39). Utilizing 1000 mb as 

-4 -1 

the pressure scale and 10 sec as the velocity divergence 
scale, these two equations can be nondimensionalized, and 
the dimensionless values of and K 2 are specified to be 
0.18. 


Following Kuo, the temperature of the cloud is regarded 
as the temperature of the moist adiabat through the condensa- 
tion level representative of the surface layer. Calculations 
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of Q and M are shown in Figs. 38 and 40. Also, for the 
c c 

modified scheme, the precipitation rate is 0.0822 in hr 
the cloud tops are at 220 mb, and the maximum of the latent 
heat released is at 400 mb. Compared to the observations, 
the deviations are only a few percent. 

The latent heat released by Kuo's original scheme is 
shown in Fig. 38. It is about the same amount as accounted 
for by the first component of the modified formulation. As 
mentioned by Ceselski (1973) the maximum latent heat released 
computed from Kuo's scheme, is often in the lower troposphere 
this is verified in Fig. 38. 



CHAPTER VI 


SUMMARY AND RECOMMENDATIONS 

The most important feature of the variational opti- 
mization scheme developed in this study is that the unknown 
upper boundary values of the o) field become solvable by 
consideration of the spatial distribution of upward and down- 
ward motions. Also, because the functional includes the 
velocity field and the u field, the procedure insures that 
both fields are optimized. Each condition and procedure is 
physically understandable and mathematically rigorous; there- 
fore, the scheme can be assumed to be realistic with a high 
degree of confidence. 

In this study, AVE II data are utilized to verify the 
reality of the optimization scheme. However, other data sets 
should also be tested. In order to provide dynamically con- 
sistent fields of atmospheric variables as initial conditions 
in a numerical forecast model utilized to verify the applic- 
ability of the modified parameterization scheme, other 
dynamical and kinematical conditions should be incorporated 
into the optimization formulation. For example, appropriate 
advection equations for dry static energy and moisture might 
be considered. 

The most important improvement of the modified para- 


40 



41 


meterization scheme is that it models physically realistic 
processes which may be verified by examination of the large- 
scale heat and moisture budget equations. Another significant 
feature is that the modified scheme is consistent with a 
variable cloud ensemble model which may be time dependent or 
static, and may or may not account for entrainment, detrain- 
ment, and/or evaporation. If large-scale moisture convergence 
above the cloud base is neglected and cloud mass flux is 
assumed to be constant with height, the latent heat released 
in the modified scheme becomes proportional to the cloud mass 
flux at the cloud base and is similar to that proposed by 
Ooyama (1969) in his numerical simulation of the life cycle 
of tropical cyclones. Furthermore, the vertical distribution 
of the latent heat released, computed from the modified scheme, 
has a maximum value at the 400-mb level which is physically 
understandable and is verified by observations (Reed and 
Recker, 1971). These features are not present in Kuo's 
scheme. 

Although results obtained with the modified scheme 
are very encouraging, they are still far from perfect. For 
example, the two-layer model which is utilized to determine 
the cloud mass flux should be improved. Extensions of this 
study probably should be focused on the following subjects. 

(a) Improved treatment of evaporative downdraft 

cooling . This may be the biggest defect of current 
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versions of all parameterization models. The down- 
draft cooling is a very effective mechanism to trans- 
port the horizontal momentum of upper levels downward 
and increases the low-level mass convergence around 
the cloudy area. The low-level mass convergence 
induced by downdraft cooling can not be observed by 
the present observational network and should be para- 
meterized in terms of large-scale variables and bulk 
properties of a cloud ensemble. This may be performed 
by dividing the cloudy area into downward and upward 
regions (Asai and Kasahara, 1967) ; the interaction 
between these two regions and between cloudy regions 
and the environment may be revealed by consideration 
of two sets of cloud ensemble equations. Combination 
of considerations of downdraft cooling and of the 
spectral distribution of cloud properties is a very 
challenging problem for future research. 

(b) Improved treatment of the sub-cloud layer . It 
is evident that accurate determinations of the mass 
flux, temperature, and moisture content at cloud base 
are very necessary for accurate computations of latent 
heat released and of bulk properties of the cloud 
ensemble. However, because of the poor understanding 
of atmospheric turbulence, estimates of the flux of 
heat, momeritiam, mass, and moisture into the cloud 
base from the planetary boundary layer are still 



highly uncertain. Variational optimization of ther- 
modynamic properties and of the wind field may be 
utilized to improve the realiability of the computa- 
tions by employing a planetary boundary layer model. 

A feasible model should be consistent with the large- 
scale motion field above the planetary boundary layer 
such that the derived eddy heat flux, Y, is nearly 
zero at the top of atmosphere (or at the tops of the 
clouds) . 

( c ) Inclusion of the spectral distribution of cloud 
properties . The clouds may be divided into continu- 
ous or discrete categories according to the heights 
of the cloud tops ; it may be assumed that the cloud 
bases are the same. Since the cloud properties are 
different for different categories, the latent heat 
released should be different for different clouds. 

In order to find the cloud properties and the latent 
heat released, the spectra of the cloud population 
should be found first. This may be the key problem 
for further research. 

(d) Improved treatment of radiation term . The radia- 
tion term in this study is computed from Rodgers' 
climatological data. Indeed, clouds can change the 
radiation dramatically. Some tests show that clouds 
in an overcast situation change the radiational 
cooling by one order of magnitude. However, better 
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radiation results depend on better cloud information, 
and the cloud information is vinknown in the convec- 
tion-parameterization problem. In order to solve 
this problem, combination of the parameterization of 
cumulus convection and a radiation model may be 
necessary. 

(e) Incorporation of scheme into a numerical pre- 
diction model . After applying the modified paramet- 
erization scheme in a diagnostic sense on several 
data sets, necessary modifications should be deter- 
mined in order that the scheme may be incorporated 
into a regional numerical forecast model such as 
Kreitzberg's (1974). The scheme should then be 
inserted into the forecast model and tested on real 


cases'. 
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APPENDIX A 


KUO'S PARAMETERIZATION SCHEME 


The large-scale flow variables are interpolated to the 

grid points of a horizontal grid system with grid length 

Ax = Ay = A. The value of the large-scale flow variable x^j 

at any grid point can then be taken as representing the 

2 

average of x over the area A = A centered at this point, 
while the actual value of x is given by the Siam of the aver- 
age value and the departure, i.e., 

X=X+X'. x=r/xdA, x'=0. (A.l) 

A •'A 

The grid area A is chosen to be much larger than the area 
occupied by a single cumulus cloud and its surrounding des- 
cending region such that a large number of clouds are included 
in A. 

The equations for the potential temperature 0, the water 
vapor mixing ratio q, and the horizontal velocity ^ of the 
large-scale system can be written in (x, y, p, t) coordinates 
as shown below. The equation for energy conservation is 


de 

dt 



Ln p _ Lu ^ 

C C ^1 

P P 



V*^'0' , 


(A. 2) 


Also, the equations expressing conservation of moisture and 
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momentum are 

K S - ’'q = -=1 - If 

^+ficx^+V$-? = -V^' , (A. 4) 

at op 

respectively, where and are the condensation rates pro- 
duced by the large-scale motions and by the subgrid-scale 
convective motions, respectively, L is the latent heat of 
condensation, is the heating rate by radiation and turbu- 
lent diffusion, and ? are the rates of turbulence diffusion 
of moisture and momentum, respectively, d/dt is 3/3 t + V • V 
+ aj|— , and tt = (P/p)^'^^p. The other symbols have their usual 

dp 

meanings . 

Assuming that at any given moment the active clouds , 
including their strongly descending parts, occupy the frac- 
tional area a, while their environments occupy the fractional 
area (1-a) , and denoting the flow variables in the active 
cloud regions by a subscript c and those in the surrounding 
regions by a subscript d, we then have 

X = (l-a)x^j + ax^ / (A. 5) 

where x stands for either w, V, 6, or q. Also, by means of 
algebraic manipulations we obtain 

t^'X' = 2 ^ - w) (x^, - X) . (A. 6) 

Since a « 1, to >> w, we get 

c 
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w'x' * - > 

- ^c^Xd-Xc> 'if Xd = X . 


(A. 7) 


Deep ciimulus towers and cumulonimbi appear to be associated 
with a deep conditionally unstable layer and the presence of 
large-scale convergence (Riehl 1950, Riehl and Malkus 1961). 
The former of these two conditions make it possible for huge 
cumuli to penetrate into the upper troposphere and the lower 
stratosphere, while the latter condition provides a general 
lifting mechanism to trigger the convective instability. 
These two conditions can be represented as follows; 


HiH2A9e > 

" > ^2^Ps ■ P*> 


(A. 8) 


where A9^ is the maximum difference of the equivalent potential 
temperature in the conditionally unstable layer, is the 
depth of this layer, H 2 the height difference between the level 
where 9 is a minimum and the level above when 9 is equal to 
its maximum value below, t is the period of the large-scale 
flow, p -p* is the lift needed for the surface air at the 
level Pg to become saturated, the maximum low-level p- 
velocity, and and K 2 are critical values of instability 
and low-level convergence, respectively, which ensure that 
deep convergence will develop. Kuo used the net convergence 
of moisture into the vertical column of air of unit cross 
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section produced by the large-scale flow and by evaporation 
from the ground as one f undcimental parameter. Let this quan- 
tity be Mj^; then. 


= , ,.. 9 , 

where subscripts s and o indicate values at the surface and 
at a nearby higher level, respectively, and Cj^ is the drag 
coefficient. 

The cimount of moisture needed to create a deep cumulus 

of area a and pressure depth with temperature and 

saturation mixing ratio q (T ) is 

o o 


aMj^ 



T) + (q^ - q)] dp . 


(A. 10) 


This is supplied by the large-scale convergence of moisture 
in time t. Hence, 


a 



(A. 11) 


where t is the half-life of the cloud (approximately 30 min- 
utes) . 


Let's assume that a fraction (1-b) of the total conver- 
gence of moisture is condensed and either precipitated out 
as rain or was carried away, while the remaining fraction b 
of is stored in the air to increase the humidity (including 
the influence of evaporation of condensed water) . That is to 
say , we have 



(A. 12) 
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Now, according to Kuo (1965) , the rate of the release of 

latent heat Q is 
c 

Qc = 7 Cp(T^ - T) (1 - b) (A. 13) 


Also, according to Kuo (1974), the vertical velocity u is 

o 

approximated as 



(A. 14) 


where 06^/3?)^^ = 30/3p + ^(3q /3p) . 
e er p 

Kuo suggested that the temperature of the cloud can be 
taken as the temperature T of the moist adiabat through the 
condensation level of the representative surface air, as a 
first approximation. 



APPENDIX B 


RAWINSONDE STATIONS FOR AVE II PILOT EXPERIMENT 



Fig. B.l, Rawinsonde stations for AVE II Pilot Experiment. 
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List of Rawinsonde 


Station Number 

11001 (MSF) 

22001 (OUN) 

22002 (FSI) 

22003 (LNS) 

22004 (FTC) 

22005 (CHK) 

201 (EYW) 

202 (MIA) 

208 (CHS) 

211 (TPA) 

213 (AYS) 

221 (VPS) 

226 (MGM) 

232 (BVE) 

235 (JAN) 

240 (LCH) 

248 (SHV) 

250 (BRO) 

255 (VCT) 

260 (SEP) 

261 (DRT) 

265 (MAF) 

304 (HAT) 

311 (AHN) 

317 (GSO) 

327 (BNA) 

340 (LIT) 

349 (UMN) 

363 (AMA) 

402 (WAL) 

405 (lAD) 

425 (HTS) 

429 (DAY) 

433 (SLO) 

451 (DOC) 

456 (TOP) 

486 (JFK) 

494 (CHH) 

518 (ALB) 

520 (PIT) 

528 (BUF) 

532 (PIA) 

553 (OMA) 

562 (LBF) 

606 (PWM) 

C37 (FNT) 

645 (GRB) 

654 (HUR) 

655 (STC) 

662 (RAP) 

712 (CAR) 

734 (SSM) 

747 (INL) 

764 (BIS) 


Table B.l 

Stations for AVE II Pilot Experiment 
Location 

Marshall Space Flight Center, Alabeuna 

Norman, Oklahoma 

Ft. Sill, Oklahoma 

Lindsay, Oklahoma 

St. Cobb, Oklahoma 

Chickasha, Oklahoma 

Key West, Florida 

Miami, Florida 

Charleston, South Carolina 

Tampa, Florida 

Waycross , Georgia 

Eglin AFB, Florida 

Montgomery , Alabama 

Boothville, Louisiana 

Jackson, Mississippi 

Lake Charles, Louisiana 

Shreveport, Louisiana 

Brownsville, Texas 

Victoria, Texas 

Stephenville , Texas 

Del Rio, Texas 

Midland, Texas 

Hatteras , North Carolina 

Athens, Georgia 

Greensboro, North Carolina 

Nashville, Tennessee 

Little Rock, Arkansas 

Monette, Missouri 

Amarillo, Texas 

Wallops Island, Virginia 

Dulles Airport, Virginia 

Huntington, West Virginia 

Dayton , Ohio 

Salem, Illinois 

Dodge City, Kansas 

Topeka, Kansas 

Kennedy Airport, New York 

Chatam, Massachusetts 

Albany , New York 

Pittsburg, Pennsylvania 

Buffalo, New York 

Peoria, Illinois 

Omaha, Nebraska 

North Platte, Nebraska 

Portland, Maine 

Flint, Michigan 

Green Bay, Wisconsin 

Huron, South Dakota 

St. Cloud, Minnesota 

Rapid City, South Dakota 

Caribou, Maine 

Sault Ste Marie, Michigan 

International Falls, Minnesota 

Bismarck, North Dakota 



APPENDIX C 


SOLUTIONS OF EQUATION SET 

AND THEIR RESPONSE FUNCTIONS 


As shown in Chapter III/ the fiinctional is defined as 
1 =j j + (v-v)^ + + 3^ (Vo)^* Vo) ^) ] dp dy dx. 


X y p 

and is associated with two strong constraints; 

= oj_ - f V*V dp , 


t = “S 

y'y'o)^ dx dy = 0 , 


(C.l) 

(C.2) 

(C.3) 


X y 


where each symbol is defined in Chapter III. 

Let us take the variation of the functional I and set 
it to zero , i . e. , 


61 = ij'ff {[(u-u) - 33^1 "Su + [ (v-v) - 3^] < 5 v + [a^ + 

X y p 

A /* X” X 

+ ^2 “ 6o)^}dp dy dx + J J' 

t s 


y p 

r r ^"^2 r r -2 ^“t ^"^2 

/ /‘e" 83 TIx=x" '>1' 

X p y p 


^// 

X p 


,2 3“t y=^2 


[3 


dp dx = 0 / 


3y --y^y 


(C.4) 


where x^^ and X 2 are the boundary values of x corresponding to 
the value of y and p held constant in the x-integration, y^^ and 
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are the boundary values of y corresponding to the value of 

. 2 

X and p held constant in the y-integration, and V xs a two 
dimensional Laplacian operator. Then, the associated Euler- 
Lagrange equations are: 

u = u + ^ , (C.5) 




~2 

a 0 ). + 


V = V + 

X, 


3Xj^ 

3^ ' 


t Po-Ph 


= X., - 


a 0)^ 


(C.6) 

(C.7) 


For the elimination of the last four terms in (C.4), there 
are several choices of boundary conditions associated with 
Lagrange multiplier X^, the derivatives of and the varia- 

tions of u, V, and o)^. In this study, the associated boundary 
conditions are: 

Xj^ = 0 , (C.8) 

n*7u^ = 0 , (C.9) 

where n is an unit normal vector along the boundary. The 
physical explanations of the variational formulation and the 
associated boundary conditions are presented in Chapter III. 

Eqs. (C.2), (C.3), (C.5), (C.6), and (C.7) are five 

equations for five unknowns: u, v, o)^, X^, and X 2 . Substi- 

tuting (C.5) and (C.6) into (C.2), we get 

“t “ \ + <Ps ■ • 


(C.IO) 
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The substitution of (C.IO) into (C.7) leads 

eV[(p^ - p^)7%] - a^(p^ - Pt>v2x^ + ^2 - 

s ^ 

(C.ll) 

From Gauss' theorem and the boundary condition (C.9), we have 

V^o)^ =JJ (V^o)^) dx dy/y* j ' dx dy 
X y X y 

= J (n*Vu^) ds/ J' y*dx dy 
C X y 

» .0 , 


where C is the boundary of the area considered, and ( ) is 
an area-average operator. After applying the atea- average 
operator to (C.7), we get 




(C.12) 


where the strong constraint (C.3) has been utilized. 

Eqs. (C.ll) and (C.12) contain two unknowns: and X 2 * 

Solutions can be obtained by an iterative technique. Fortu- 
nately, because the equations are linear and X 2 is constant, 
the solutions can be determined without using an iterative 
procedure. In order to clarify the procedures utilized to 
obtain solutions, let us define a linear differential operator 
L as : 

L = sV(p^ - p^)v2 - a=(p^ - p^JV^ + p4p- • 

s ^ 


Then, (C.ll) can be rewritten as 


(C.13) 
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where l”^ is the inverse operator of L. 

Setting 

Xi = Xfa + <j) , (C.14) 

where 

= L~^ (C.15) 


and the boundary values of and (p are chosen to be zero, 
then, (C.13) and (C.12) become 


<f) = L ^(X2) , 


(C.16) 


Xo = C + (— ^) / 


Ps'Pt 


(C.17) 


where 


C = 52 5^ + ( 5 ^) 


From (C.15), Xj^ can be solved by using a relaxation method. 
After 0 and X 2 are determined from (C.16) and (C.17), the 
complete solution is obtained from (C.14). Eqs. (C.15) and 
(C.16) show that Xj^ is a part of the solution of X^ which 
corresponds to the curvature of the w^-field, and (p is the 
other part of the solution of Xj^ which corresponds to the 
mean value of These properties will be discussed in 

detail later. 

In order to employ iterative procedures to solve for 
(J) and A 2 , let us rewrite (C.16) and (C.17) as 


^(k) ^ l"^(X^^^) 


f 


(C.18) 
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= c . 


<p 


(k-1) 


Ps"Pt 


(C.19) 


where (k) indicates the k-th iteration. Let zero be the first 
guess for the solution of (p , and define 


L~^(C) 

Ps'Pt 


aC , 


i. e. , 


a = 


L ^(C) 

Ps-Pt 


Then, (3.21) and (3.22) become 


= c z 

k=0 


k 

a , 


(C.20) 




(n) 



n 

E 

k=0 


k 

a , 


and 



L ^(C) . 


(C.21) 


(C.22) 


The derivation of (C.20) and (C.21) are mainly based on the facts 
that the operator L is linear and is constant. It is obvious 
that if |a|< 1, the iteration procedure is convergent provided 
C and are bounded. Since a is a function of 5 and 3/ a 

diagram, as shown in Fig. C.l, can be constructed to determine 
the convergent and divergent region in (a, 3) space. Conver- 
gence of the iterative procedure has been investigated for 
several cases. If the variables are nondimens ionali zed by 
suitable scales (e.g. the scale for the vertical p velocity 
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should be 10 mb sec , and the mass convergence scale should 

be 10 ^ sec etc.), all tested cases show that the iteration 

procedure is convergent (rapidly in most cases) for a wide 
range of values of a and 3. 

In the convergence region, the solution of (C.20) and 
(C.21) is 


A 


2 


^2 


d 

1-a ' 






1-a " 


(C.23) 


and the complete solution of (C.ll) and (C.12) is obtained by 
use of (C. 14) . 

Following the same concepts utilized in obtaining solu- 
tions to (C.14), (jii^ is written as 

0)^ = + 0)^'*’ (C.24) 

where corresponds to and u)^'^ corresponds to <J). Referring 

to (C.IO), and may be chosen as 

“t^ = Wt + (Pg - ' (C.25) 

= (Pg - , (C.26) 

and the associated boundary conditions are 

= n-Vo)^*^ = 0 . (C.27) 


Boundary conditions different from (C.27) may be posed, al- 
though we prefer (C.27) because they are simple and are 
consistent with (C.9). According to (C.25) and (C.26), 
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and can be determined from and (p, and then can be 

obtained from (C.24). 

Figs. 9-26 provide examples of this optimization 
scheme. Comparisons of these vertical velocity fields before 
and after optimization are discussed fully in Chapter V. 

In order to better tanderstand how to choose suitable 
values for a and 3, the associated response function has to 
be found. Let us assume that 

Ps ' Pt ' Po P' ' 


00 

• ' = s 


he ] 

whole region, and p' is the fluctuation of (p^ - p^) from p^. 


k=l 

where e is a small parameter, p^ is the mean value over the 


Since the ratio of p ' to p is smaller than 0.1 and is less 

o 

than 0.02 at most grid points, p' is treated as a small per- 
turbation. The differential operator L can be expanded as: 

L = (3^p - a^p + — ) +{3^V^p'V^ - a^p'V^ 

o o Pq 


— E (- = 

Po k=0 Po 


+ 0(e) , 


where 


^o = ^^Po^"^ - ' 


and 0(e) indicates the order of magnitude of e(less than 0.1 
in this study). To the zeroth order approximation, (C.15) 
and (C.16) become 



The solution of 4» is simply 


<|» = PqX2 , (C.29) 

From Fourier transformation (FT) theory, any function 
X(x,y) and its Fourier transform C(m,n) are represented as 


oo oo 


X(x,y) = FT~^(C(m,n)) = Z Z C (m, n) e^ ^ dm dn , 

n=-oo m=-“ 


L 

y X 


C (m,n) = ft (X ( x 




-i (mx+ny ) 


■_x -_Ji 
2 2 


(C.30) 

dx dy , 
(C.31) 


where m, L^, n and are the wave numbers and the width of 
the domain along the direction of x and y, respectively. 
Utilizing (C.30) and (C.31), (C.28) becomes 

^b ^ C^ (m,n) 

C^ (m,n) 


(6^k^+a^)k^p +— 

° Po 

2 2 2 
where k = m + n , 

C^(m,n) is the Fourier transform of Xj^, and C^(m,n) is the 
Fourier transform of Similarly, from (C.25), we have 
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C^(m,n) 


= 1 - 


p 


(B^k^+a^)k^P^+^ 


where C^(m,n) is the Fourier transform of and C^(m,n) is 

the Fourier transform of o)^. 

b b 

In spite of the value of a, approaches 1 and 
approaches zero as B approaches zero for any m, and n. It 
implies that c“(m,n) approaches C^(m,n) and Aj^ approaches zero, 
i.e. , no filtering of the field. It is always true even if 
p' is not negligibly small. On the other hand, as B becomes 
very large, we have 

b 1 
Y ~ 


and 


Po^' 


~ 0 
'w 


The above equations show that the values of Yj^ ^.re dominated 
by long waves as B becomes large. It also implies that 
approaches zero as B approaches infinity. 


Since X 2 is a constant, the combination of (C.29) and 


(C.26) implies 


wj = 0 , 


identically. The above equation indicates that when p' is 
small, most of the adjustment of the vertical velocity (or mass 
convergence) is due to the part. It will be shown, in 


Ill I Hill llll■lllll II iini 


iiiiiiii III 


II I II 


■I ■ I ■ II 
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Appendix D, that when 6 is zero, the adjustment is very 
small, provided a reasonable value of a (e.g. , a = 0.1) is 
utilized. 

According to the above discussion, several numerical 
examples of the response function are shown in Fig. 42. 
Generally speaking, for synoptic-scale systems, errors in 
the velocity field are of the order of 10 percent and errors 
in the u field may be ten times larger. In other words, the 
reliability of the velocity field may be ten times larger 
than that of the w field, and, therefore, a may be chosen as 
0.1 in the optimization scheme. Since the wavelength of 
resolvable waves is about 600 km and the width of the area 
considered is about 2000 km, waves whose wave number is less 
then three should be suppressed. Also, for the synoptic 
situation under study, the distance between adjacent extrema 

I 

in the vertical motion field is between 700 and 1000 km. 
Therefore, the most important waves are those of wave number 
one or two. Values of 6 between 0.1 and 0.3 should be very 
good in this situation. 



APPENDIX D 



SOME SPECIAL CASES OF THE VARIATIONAL 
OPTIMIZATION SCHEME 


If 3 is zero, (C.17) becomes 

~ 2 


= 0 


(D.l) 


“ (Ps-Pt^ 


~2 2 

Since a and (p - p ) are both positive, the solution of 

S u 

(D.l) may be expressed in terms of the Green's function, 
G(x,y,C,n) as 

Xj^(x,y) = y* (x,y ,c,n) ' 

where n and s are the coordinates perpendicular and tangential 
to the boundary C, respectively. Because is zero along the 
boundary, the solution of (D.l) is simply 

Xj3 = 0 , (D.2) 

identically. As a approaches zero, (C.18) and (C.IO) become 

Xi = <|) = X2(Pg-p^) (D.3) 

and 

Substituting (D.4) into (C.3), we get 


X 


2 



I(P3-Pt>V^(Ps-Pt) 


(D.5) 
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(D.3), (D.4) and (D.5) imply that if (p^ -* p^) is a con- 
stant, ( 0 ^ is equal to In this case, the adjustment of 

5^ is impossible. In other words, the horizontal curvature 
of the (p - p ) field is the only mechanism to adjust the ob- 

S o 

servation field such that the continuity equation and the global 
upper boundary condition of vertical p-velocity are satisfied. 

If 3 is zero and a approaches infinity, is zero 
identically, and (C.18) and (C.IO) may be written as 


2 2 
V'^X, = - ^ 

1 (Pg-Pt) 

“t “ “t ■ ^2 • 


(D.6) 


The above two equations imply 


Since 


and 


we have 


^2 = “t 


(D.7) 


0) = (1). 


s 

0) - 3 = - / (V*V-V*V) dp 

•^Ps 


-I 


- I dp , 

® JPs 

P 


V*V dp , 


= V^X^ dp 


Ps"P ~ 

= - “t 


(D.8) 


s *■ t 

(D.8) is the differential form of O'Brien’s formulation (1970) 
Because the correction velocity field is irrotational , the 
optimized wind field can be determined after the correction 
velocity potential is determined from (D.6) and (D.7). 
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If 6 approaches, infinity , then for any boianded a, (C.15) 
implies 

= 0 , 

provided and are bounded. Siibstitution of the boundairy 
condition (C.27) leads to 


= constant. 


Combination of (D.9) and (C.25) yields 
y2^ _ _ ^t constant 


Ps'Pt 


Ps‘Pt 


(D.9) 


(D.IO) 


Furthermore, from (C.16) we have 

V^[(Pg-p^)V^(t.^°^] = 0 , 


(D.ll) 


provided Aj^ is bounded. The substitution of (C.23) into 
(D.ll) leads to 

V^[(Pg-p^)v2({,] = 0 , 

or 

V^a)^*** = 0 , 

by the use of (C.26). From the boundary condition (C.27), we 
have 

( 0 ^*^ = constant , (D.12) 


or 


y2. _ constant 

Ps~Pt 


where d and a are bounded provided a is bounded. Substitution 
of (D.9) and (D.12) into (C.24) leads to = 


constant. Then, 
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applying the boundary condition (C.9), we get = 0 , iden- 
tically. Therefore, Eq. (C.IO) reduce to 




P0-P4 


The above equation is exactly McGinley's formulation (1973). 




Fig. 1. Data grid and interior working grid 


















Fig. 7. Same as Fig. 3 except for 300 mb 
with contours drawn at 120 m intervals and 
isotherms constructed at 4°C intervals. 
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( 230 ) — — 

Fig. 8. Same as Fig. 3 except for 200 mb 
with contours drawn at 120 m intervals and 
isotherms constructed at 4®C intervals. 


Fig. 9. Vertical velocity, u = dp/dt, 
at surface before optimization. Iso- 
lines are drawn at intervals of 0.4 yb 
sec"^. Positive values refer to decent 
negative values refer to asceht. 
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Fig. 10. Same as Fig. 9 except for 
800 mb with isolines drawn at intervals 
of 2 yb sec~l. 



Fig. 11. Same as Fig. 9 except for 
700 mb with isolines drawn at intervals 
of 2 yb sec"!. 
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Fig. 12. Same as Fig. 9 except for 
600 mb with isolines drawn at intervals 
of 4 yb sec"”^. 



Fig. 13. Same as Fig. 12 except for 
500 mb. 




Fig. 14. Same as Fig. 12 except for 
400 mb. 



Fig. 15. Same as Fig. 12 except for 
300 mb. 








Fig. 16. Same as Fig. 12 except for 
200 mb. 



Fig. 17. Same as Fig. 12 except for 
100 mb. 
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Fig. 18. Optimized vertical velocity, 
omega, at 800 mb with isolines drawn 
at intervals of 0.2 yb sec“l. 



Fig. 19. Same as Fig. 18 except for 
700 mb with isolines drawn at intervals 
of 2 yb sec”l. 





Fig. 20. Scune as Fig. 19 except for 
600 mb. 



Fig. 21. Same as Fig. 19 except for 
500 mb. 
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Fig. 22. Same as Fig. 19 except for 
400 mb. 



Fig. 23. Same as Fig. 19 except for 
300 mb. 





Fig. 24. Same as Fig. 19 except for 
200 inb. 



Fig. 25. Same as Fig. 18 except for 
100 mb with isolines drawn at intervals 
of 0.02 yb sec“^. 





Fig. 26. Horizontal divergence at 
test point. 



Fig. 27. Vertical velocity/ w = 
at test point. 
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Fig. 28. Observed temperature and dew-point tempera- 
ture at test point. 







Pressure (mb) 


O„(^C(fay-0 


/li 

h^'h 

// 

// 

( » 

\ 

• \ 

w 

\ \ 

■\\ 

‘•••J I 

80 85 

Cah gm“^ 

Fig. 29. The dry static energy, s 
(solid) , moist static energy, fi 
(dotted) , and saturation moist static 
energy, fi* (dashed), of environment. 
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Fig. 30. The apparent heat source, 
Q, (solid) , the apparent moist sink, 
Q 2 (dashed) , and the radiational 
heating from Rodger's results (1967) 
(dotted) . 
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Fig. 31. The derived vertical eddy 
heat flux, Y. 
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Fig. 33. Released latent heat, com- 
puted from large-scale budgets and 
Kuo's formulation. 
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Fig. 34. Cumulative mass flux distri- 
bution within clouds, (solid), and 
the corresponding environmental mass 
flux, M (dashed) , computed from the 
large-scale budgets and Kuo's formula- 
tion of the latent heat released by 
subgrid-scale convection. 
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Fig. 36. The excess mixing ratio, 

computed from the large-scale 
budgets and Kuo's formulation of the 
latent heat released by subgrid-scale 
convection. 
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Fig. 37. Cumulative mass flux distri- 
bution within clouds, (solid), and 
the corresponding environmental mass 
flux, M (dashed) , computed from the 
large-scale budgets and the modified 
formulation of the latent heat re- 
leased by subgrid-scale convection. 



Fig. 38. The latent heat released 
by subgrid-scale convection, com- 
puted from (a) Kuo's scheme (dotted), 

(b) the modified scheme (dashed) , and 

(c) the large-scale budgets and the 
modified formulation (solid) . 
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Fig. 39. The excess temperature, 
- T (solid) , and the excess 
mixing ratio, ^ (dashed) , 

computed from the large-scale 
budgets and the modified 
formulation of latent heat 
released. 
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Fig. 40. Cumulative mass flux distri- 
bution within clouds, (solid) , and 
the corresponding environmental mass 
flux, fl (dashed) , computed from the 
modified scheme. 
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Fig. 41. The first (solid), the second 
(dashed) , and the third (dotted) compon- 
ents of the expression for the latent heat 
released, computed from modified scheme. 
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Fig. 42. The response function for 
a = 0.1. 
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